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ABSTRACT 

We investigate the structure of accretion disks around massive protostar applying steady 
state models of thin disks. The thin disk equations are solved with proper opacities for dust 
and gas taking into account the huge temperature variation along the disk. We explore a wide 
parameter range concerning stellar mass, accretion rate, and viscosity parameter a. The most 
essential finding is a very high temperature of the inner disk. For e.g. a 10 Mq protostar and 
an accretion rate of ^ lO~'*M0yr~^, the disk midplane temperature may reach almost 10^ K. 
The disk luminosity in this case is about 10'* and, thus, potentially higher than that of a 
massive protostar. We motivate our disk model with similarly hot disks around compact stars. 
We calculate a dust sublimation radius by turbulent disk self-heating of more than 10 AU, a 
radius, which is 3 times larger than caused by stellar irradiation. We discuss implications of this 
result on the flashlight effect and the consequences for the radiation pressure of the central star. 
In difference to disks around low mass protostars our models suggest rather high values for the 
disk turbulence parameter a < 1. However, disk stability to fragmentation due to thermal effects 
and gravitational instability would require a lower a value. For a = 0.1 we find stable disks 
out to 80 AU. Essentially, our model allows to compare the outer disk to some of the observed 
massive protostellar disk sources, and from that, extrapolate on the disk structure close to the 
star which is yet impossible to observe. 
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Massive stars play a vital role in order to un- 
derstand the dynamical evolution of clusters in 
which they are the major source of heavy ele- 
ments and UV radiation. During their short life 
time, they impact their surrounding by a num- 
ber of physical processes such as jet like out- 
flows, strong winds, photo evaporation, expand- 
ing HI I regions and eventually s upernova explo - 



(|Beuther et all l2002t iBeuther fc ShepherdI [2005h . 
although the launching mechanism for outflows 
have been mainly investi gated for low mass pro- 



tostars so far (see e.g. Casse fc Keppend 2002 : 



Pudritz et al.ll2007t lFendt]l2009[ ). 



sions (jZinnecker fc Yorkd 120071 : iBeutheil 120071 ). 

Collimated molecular outflows and jets from 
young massive protostars have been observed 



''Member of the International Max Planck Research 
School for Astronomy and Cosmic Physics at University 
of Heidelberg, IMPRS-HD 



Understanding massive star formation has been 
a very active field of research for both the ob- 
servers and the theorists. The basic differences of 
high and the low mass star formation process are 
that of timescales, energy scales, mass flow rates 
and high luminosity of the central star. Low mass 
stars have well defined phases in their formation 
process and they only start burning h ydrogen af- 
ter ac cretion of all matter is done (Sta hler fc Palla 
20051 ). In contrast, massive stars have very short 



Kelvin Helmholtz time of ~ 10 — 10 yr, and thus 
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start burning hydrogen even when they are still 
accreting. One of the problems in the formation 
of massive stars is that, at some point in time, the 
strong radiation pressure from the luminous cen- 
tral star may exceed the Eddington limit in spher- 
ical symmetry and does not allow the matter to 
accrete anymore. This seems to limit the mass of 
the star to be fo rmed via the spherical accretion 
to about 20 Mq (jWolfire fc Cas sincUi 1987). 

However, disk accretion may add dynamical 
pressure to the accreting matter which may help 
overcoming the stellar radiation pressure. Thus, 
the star formation scenario for low and high mass 
stars could in principle be the same (e.g. Ketcl 
20071 ) ■ suggesting that high mass star formation as 
a scaled up version of the low mass star formation 
process. With our paper we aim to construct a 
global model (length scales of ~ O.lAU to lOOAU) 
of the accretion disk around massive young stars. 
By fitting the outer disk structure to the obser- 
vations, this will allow us to investigate, in follow 
up projects, the feasibihty of certain physical pro- 
cesses in inner disk which could be essential for the 
pre-stellar evolution as angular momentum trans- 
port or processes responsible for launching of out- 
flows. 

Our paper relies on one essential assumption 
- that is the disk accretion rate. Both estimates 
following the observed outflow rates and the time 
scale of mass aggregation suggest accretion rates 
in the range of 10" ^-10~'^Mf7^vr~^ (iBeuther et al 



2002t IZhandl2005t ICrave &: Kumaij|2009t ). How- 
ever, one of the aims of our study is to understand 
whether such values are compatible with theoret- 
ical accretion disk models. 

In the following section of this paper, we discuss 
about the modelling of the disk structure using the 
standard disk equations, in §3 we describe the disk 
opacity model which has contribution from both, 
dust and gas and in §4 we present our results and 
give estimates of the sublimation radius and then 
in §5 discuss about the stability of the disk and 
how the viscosity parameter a has an effect on 
it. Also in §6 we compare our results with the 
recent observations of disks in massive young pro- 
tostars. Finally, in §7 we present our conclusions 
from this simple approach to study disks around 
massive stars. 



Accretion Disk Model 



In order to study the global disk structure 



(Shakura & Svunvaev 


1973: MorfiU & Woodlll989: 


Ruden & Pollackll991 


; Frank et al.lll992; Steoinski et al. 


19931: Bell & Lin 199^ 


1) with appropriate modifi- 



cations. The essential part of our model is to take 
into account the proper gas and dust opacities for 
a huge temperature regime as indicated for high 
accretion rates on to mass ive stars. The origi- 



nal st andard disk model bv I Shakura fc Svunvaev 
(jl973[ ) relates to solving the stationary state hy- 
drodynamic equations for the thin disk assuming 
certain opacity laws for different regions in the 
disk. The disk is classified to have three regions 
(A,B, and C) depending upon the opacity and 
the pressure. The region C is the outermost re- 
gion where the gas pressure is much greater than 
the radiation pressure in the disk and free-free 
opacity dominates in this region, whereas region 
B is the one in which gas pressure is still greater 
but the Thomson scattering is now dominating 
over Kramer's opacity. The region A is the one 
which is closest to the central object and here 
the radiation pressure in the disk is much greater 
than the gas pressure. The thin disk model was 
extended to include dust opacities and has been 
applied t o disks around low mass stars and solar 
nebula. (IRuden fc Pollackl Il99ll: IStepinski et al 



19931 : IStepinskil Il998l : iDel Popolo fc Eksil l2002h . 

The application of such a thin disk model to mas- 
sive stars is a simple approach to understand the 
dynamics of the disk especially in the inner most 
region which can not be resolved with the help of 
present day telescopes. We look for steady state 
solution of the thin disk model using the hy drody- 



namic equations of the thin disk (e.g. Fra nk et al 



19921 ) . The vertically averaged surface density S is 
related to the mass density p and the scale height 
H. 

E = 2Hp (1) 

The scale height is estimated using the approxi- 
mation for thin disk 



(2) 



where fl is the Keplerian angular velocity and Cg, 
the sound speed is defined in terms of temperature 
and other physical constants like Boltzmann con- 
stant fc, proton mass nip and the mean molecular 
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weight ^. 



( k 



3. Disk opacities 



\fimp 
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(3) 



The energy balance equation which imphes that 
the amount of energy radiated from the disk is 
equal to the amount of energy produced by viscous 
heating. 



16cr 

3Ek 



(4) 



The viscosity v is defined by introducing the pa- 
rameter a. 

V — aCgH (5) 

The Rosseland mean opacity k is related to the 
mass density and the temperature in form of power 
law with kq as constant. 



kqP T 



(6) 



The condition for the conservation of angular mo- 
mentum in steady state can also be derived. 



M 



r 



(7) 



The dimensionless viscosity factor a is intro- 
duced basically to parameterize viscosity by as- 
suming that is provided by turbulent motions 
present in the disk. The physical significance of 
equation (O is that the largest size of turbulent 
eddy that can be sustained in the disk is of size 
equal to the disk scale height and maximum speed 
of the turbulent motion is that of sound speed, 
this implies values of a < 1. In order to solve 
the above equations consistently, one needs the 
form of opacity. For the opacity of the form of 
a power law given by equation ^ , one gets power 
law solutions of all the dynamical quantities in 
the disk. For instance, in case of the form of 
Kr amer opacity law (k = 6.6 x 10~^^p^r~2 ) used 
bv IShakura &: Svunvae-J (|l973l ). the central tem- 
perature in the disk around a typical massive star 
in the region where Pgas > -Prad and as > at is 



Tc(r) = 1.74 X 10'' K (a) 
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where as is the free free opacity and at is the 
opacity due to electron scattering. 



The work done bv lShakura k SvunvaevI (|l973l ) 
was mainly focused on hot accretion disks around 
black holes, the contribution to the opacity in 
their case was mainly from the electron scatter- 
ing or free-free emission depending upon the op- 
tical depth in the disk. However in disks around 
massive young stars, contribution to opacity also 
comes from dust present beyond the sublimation 
radius, many molecular and atomic lines and other 
scattering processes. Since massive stars have 
large radiation fields that heavily affect the huge 
amount of dust present in the envelope and also 
the gas that is present in the inner part, one has to 
take into account a proper contribution of opaci- 
ties from dust and gas to have a consistent accre- 
tion disk model. 

In our model the dust opacity does not have any 
explicit dependence on frequency, so the effective 
value of the opacity 



Kcff 



(9) 



simplifies to Kcff = Hu 

In our work, we follow that the matter is ac- 
creted via a disk investigating the possible loca- 
tion of the dust sublimation radius. We expect 
that the active disks will be dominating in destroy- 
ing the dust in the disk at a much further radius 
then it would have been destroyed by the heating 
from star. This would result in lowering the radi- 
ation pressure on the dust and allowing the infall 
of matter on the central star. 

Apart from some metal silicates all the dust 
grains that are present sublimates at around 
1500 K. Since at this temperature, the gas and 
the dust opacities vary substantially, it is essential 
for calculating radiation pressure to consider the 
proper opacity. 

3.1. Dust opacities 

There are many factors one has to take into 
account for getting a consistent model for dust 
opacity such as the size of the grains, distribution 
of grains in the disk and the coagulation of dust 
grains. Also the temperature of the disk varies 
over a large range and this would clearly affect 
the composition of the dust and alter the gas to 
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dust ratio typically ~ 100 in the disk. This ra- 
dial variation in the gas to dust ratio will not 
only modify the gas pressure along the disk but 
also the other dynamical disk quantities. The 
usual assumption in many of the dust models 
is that the dust and gas are well coupled and 
also the gas to dust ratio is taken to be 100. 
In general, opacity has a typical dependen ce on 
temp e rature and densit y. (e.g. iMorfill & Wooc 



19891: iRuden fc Pollackl Il99lt Beh fc Li 
Ossenko pf fc Henningj Il994 : " iHelling etal 
Scmcnov ~t al.ll2003l ). 

In the regime which is dominated by dust 
(r ^ ii*), the opacity depends very weakly on 
density. In the present work we apply for the 
dust dominated r egion the model proposed by 
iRuden fc Pollackl (|l99lf l in which the Rosselanc. 
mean opacities have been given in terms of the 
power laws of the form given in equation ([6]), in 
different regimes that have been distinguished on 
the basis of temperature and mass density (see 
Table [1]) . Using this power law form of opacity, 
equations ([I]) - (O can be reduced to a form rep- 
resenting the viscosity as a function depending on 
surface density and the radial distance as 



(10) 



This form of viscosity can be then substituted in 
equation ([7]) so that all the dynamical quantities 
namely the surface density, mass density, central 
temperature, scale height of the disk can be ex- 
pressed as power laws of mass, accretion rate, ra- 
dial distance and the viscosity parameter a. For 
the parameters of a typical massive star these 
quantities (Q) can be written 



g = c, 




(11) 



e.g. (iRuden fc Pollacklll99ll :IStepinski et al.l[l995 
Del Popolo fc Eksill2002t ). The index i denotes the 
opacity regime as given in the Table [TJ {i = 
1, 2...6), Ci are constants and Ri is chosen as a typ- 
ical radius for that regime. The dust model taken 
into account does not consider explicit wavelength 
dependence of opacity. 



3.2. Gas opacities 

The opacity of gas and dust have been esti- 
mated consistently. These opacities are usually 
listed in form of a table in two parameters space, 
temperature (T) and parameter (R) which is de- 
pendent on density in a following manner 



R 



P 



in 



(12) 



20051 ) 



( Iglesias fc Rogersl Il996l : iFerguson et al 
where p is the density and Tq is the tempera- 
ture normalized to 10^ K. These opacities also 
show some dependence on metallicity as well. 

In a regime where dust begins to sublimate and 
molecules start to form, there is a sharp decrease 
in the opacity. In this regime, dependence of opac- 
ity on density is strong, unlike the opacity due to 
dust, and in general it is very difficult to get the 
true value of the opacity in this regi on so usu 



ally linear interpolation is used (Sc menov et al 
20031 ). Due to steep turn over of the opacity 



gradient with the temperature, numerical difficul 
ties arise when the dynamical quantities are esti 
mated for this region. In the prese nt work, we take 
i nto a ccount the opacity tabl^ by iFerguson et al 



(|2005[ ). which has the temperature range of our 
relevance (500 K - lO'^'^K). For the present pur- 
pose we choose a table with hydrogen fraction X 
= 0.94 and metallicity Z — 0.06. In fact, we also in 
the inner most region obtain higher temperatures 
and so we apply in continuation the OPAL opac- 
i ty tablejl for hig her temperatures (> lO**'^*^ K) 
( Iglesias fc Rogersl Il996 ). We create equispace 
grid of two parameters log(T) and log(R) as de- 
fined above and put the values from the table in 
the grid and fit a 2D cubic spline on it to get the in- 
terpolated values of opacity for any given temper- 
ature and density, but as the regime where the gas 
and dust coexists is difficult to model numerically 
as the opacity do not converge to a unique value, 
we make a linear approximation as zeroth order 
approach and compare with the analy t ical models 
prese nt, (e.g. iRuden fc Pollackl [T99ll : iBell fcThi 
I1994)). 



t http: / / webs.wichita.edu/physics / opacity 

thttp:/ /www-phys.llnhgov/Research/OPAL/opal.html 
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Table 1: Different regimes of the opacity taken from iRuden &: PoUackl (|l99ir ) with the conditions for the 
change over and description of the main components of dust in each regime 



Regime Ki Condition Description 

1 (2 X 10"*') p°T^ T < 150K dominated by Water and Ice grains 

2 (1.15 X 10^**) p°T^** T < 180 K Sublimation of Ice grains 

3 (2.13 X 10"^) p'^T"^^^ T < 1380K (^^y'''^" mainly consists of Iron and Silicate grains 

4 (1.57 X lO"") p3/**T"i** T < 1890 K {^J^Y''^^ Sublimation of refractory grains 

5 (1.6 X 10"^) p°T° T <2620k(^j^Y'^'' Molecular and atomic lines contribute 

6 (2 X 10^") p^^^T"'' T< 3200K Molecular and atomic lines contribute 



4. Results and Discussion 

In the following we present results of applying 
above standard models of disk accretion and opac- 
ities to massive star formation. We first discuss 
the outer dust dominated disk and then the inner 
gaseous disk. 

4.1. Outer disk structure : Radial Profiles 
of Dynamical Quantities 

The main contribution of opacity in the outer 
disk is from dust. The numerical values for the 
dynamical quantities expressed as given by equa- 
tion pTjl in different opacity regimes are listed in 
Table [2] The radial profile of the disk midplane 
temperature, surface density, mass density and 
scale height for a high mass star is shown in Fig.[l] 
for a IOM0 star with accretion rate of the order of 
4.2 X 10"'' M0yr~'^ and the a viscosity parameter 
a — 1. With these parameters, the central tem- 
perature reaches 1500 K around 12 AU due to vis- 
cous heating, which causes the dust to sublimate. 
Figure [T] shows two curves for each quantity. 
The solid curve is obtained by implementing the 
opacit y power laws as given by iRuden fc Pollack 
(|l99l[ ). whereas the dashed line is obtaine d by us- 
ing the opacities bv lStepinski et aL ( 19931 ) and ex- 
trapolating it to higher mass stars. 

There are few kinks seen in the plot, which 
are due to the fact that the dust opacity model 
used comprises of different regimes and in each of 
these regimes have a different form of power law 
(see Table [1]) . These reg imes are connected using 
the standard procedure (jLin k, Papaloizoul 119851 ) 
of equating the opacity in two consecutive regimes, 



Since formation of massive stars involves large 
accretion rates the surface density and the mass 
density in the disk is higher than the lower mass 
counterparts, as more matter will be injected in 
the inner region. We find the scale height ratio 
obtained from the present model is approximately 
constant in the outer disk, with a radial depen- 



dence cx 



^0.28 



and consistent with the thin disk 



approximation. 

A least square linear fit of dynamical quanti- 
ties gives for the temperature profile a power for 
the radial distance as -0.45 and for density as - 
2.35. The surface density has a rather flat profile 
as the best fit gives the dependence of the form 
These profiles are simila r to those ob- 



tained by iRuden fc PoUackl ( 199l[ ) but here ap- 
plied for a high mass star with high accretion rate. 

The parameter values discussed above were cho- 
sen i) to have most part of the inner 100 AU disk 
to be gravitationally stable and also ii) to ensure 
that the dust in the disk sublimates at a distance 
substantially further than caused by sublimation 
from stellar irradiation, (see below) 

4.2. Sublimation of dust 

Most of the dust grains sublimate when the 
temperature in the disk reaches the critical value 
of 1500 K. The radius of the disk at which this 
value of critical temperature is reached is the dust 
sublimation radius. The disk can be heated by var- 
ious processes such as viscous heat ing, stellar irra- 



diation, convection, cosmic rays (jD'Alessio et al 
J998). However, in the present work we estimate 
the dust sublimation radius via two main processes 
namely viscous heating and stellar irradiation (see 
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Table [31). 



4-.2.1. Disk Self-Sublimation of dust 



Table 2: The various dynamical quantities 
for the massive star using the interpolated 
Rosseland mean dust opacity model given by 
Ruden fc Pollack ( 199ll ).(see equation (fTTj) and 
Table HD 







Regime 


1 Ri = 400AU 




Q 


T[K] 


Height[cYa] 


p[gem"^] 


E[gem-2] 


c. 


126.3 


9.9 X lO" 


1.3 X lO"^-'' 


258.8 








Regime 


2 R2 = 250AU 




Q 


T[K] 


Height[cm] 


p[gem"^] 


S[gem-2] 


c. 


169.8 


5.7 X lO" 


3.4 X 10"^-' 


393.1 








Regime 


3 i?3 = 50AU 




Q 


T[K] 


Height[cm] 


p[gem"^] 


S[gem-2] 


c. 


814.5 


1.12 X lO" 


4.1 X 10"^^ 


907.6 








Regime 


4 R4 = lOAU 




Q 


T[K] 


Height[cm] 


p[gcm~^] 


S[gem-2] 


c. 


1613.2 


1.4 X 10^'^ 


1.8 X 10"^° 


5126.4 








R,egimc 


5 Rs = 6AU 




Q 


T[K] 


Height[cm] 


p[gcm~^] 


S[gem-^] 


c. 


1918.3 


7.1 X 10^^ 


6.5 X 10"^° 


9282.8 




Regime 6 Re = 3AU 


Q 


T[K] 


Height[ciri 


p[gcm 3] 


S[gcm-2] 


C, 


2767.9 


3.0 X 10^^ 


3.0 X 10"^ 


1.8 X 10* 



For our massive stellar disk model, we define 
that radius where the disk temperature reaches 
1500 K as disk self-sublimation radius. Essentially, 
the disk temperature is determined by M, while 
the radial velocity is governed by a. The higher 
the accretion rate, more the gravitational poten- 
tial energy from the accreted mass is converted to 
thermal energy in the inner region, thus increasing 
the disk temperature. High mass accretion rates 
~ 10^'^ — 10"'* M0yr~^ are indirectly related to 
the formation of massive stars, indicating hotter 
disks in massive stars. 

The dust sublimation radius can be estimated 
from the temperature profile in regime 4 (see Ta- 
ble [2]) as the critical value of 1500 K is obtained 
in this regime, 



T = 1.6 X 10-^ K a 



-0.058 



M 



0.101 



lO-3]vi0yr~i 

0.080 



MY"'' / r X -0.239 
M^j (tOAu) ^''^ 



Thus the dust self-sublimation radius is 

M 



0.422 



-Rsub.disk = 11 Alia 



-0.242 



10-4]VIgyi.-l 



VlOMe 



1500K 



(14) 



4-2.2. Dust Sublimation by Stellar radiation 

The second process of disk heating is by radi- 
ation from the central star. The dust sublima- 
tion radius due to the absorption of UV radiation 
from the star by the disk is dependent on the ef- 
fective temperature and the radius of the star or 
equivalently to the luminosity of the star. We can 
estimate the relation between the temperature of 
the disk and the temperature of the star just by 
equating the flux from the star that is absorbed by 
the disk to the flux emitted by it considering it as 
black body. The disk is also assumed to be locally 
isothermal, which is an appropriate assumption for 
optically thick disks The mid plane temperature of 
the disk due to heating from central star can be 
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estimated, 



1/4 



1/2 



(15) 



(e.e. Dullemond et al. I l200lh where 6 is the an- 
gle with which the radiant flux is incident on 
the flaring disk, V's denotes the fraction ol flux 
that is absorbed by the interior and i/'i is cor- 
rection factor which accounts for the fact that 
the disk interior is not fully optically thick for 
its emission. The small correction due to back- 
warming in the disk is neglected. This correction 
main ly depends on the underlying dust proper- 
ties (jMonnier fc Millan-Gabe^l2002f) . For a dust 
sublimation temperature of ~ 1500 K, we can esti- 
mate the sublimation radius due to absorption of 
radiation from the above equation (|15p noted as 
Measure A whereas Measure B we denote as the 
sublimation radi us calculated using the standar d 
formula given bv lMonnier fc Millan-Gabe^ (|2002l ). 
In Table [3] we show the typical values of the subli- 
mation radius from these measures. For a t ypical 
B2 star of mass = lOM© and Tcff = 22000 K (jLangl 



19921 ) with ZAMS value for the stellar radius as 
6 R0, one gets a dust sublimation radius due to 
heating from stellar radiation of -Rsub,* ~ 4 AU for 
measure A and ~ 3 AU for measure B. This value 
may increase by a factor of 2 considering the bloat- 
ing star model (jHosokawa fc Omukaill200^ . where 
the star bloats up to lOOR© and the effective tem- 
perature reduces to 5000 K. 

In Table [31 we compare for different stellar mass 
with typical order of mass accretion rate, the sub- 
limation radius due to self heating from the disk 
and that due to heating from the star. This values 
are estimated for two different values of a = 0.1, 1. 
It is evident from Table [3] that the dust sublima- 
tion radius due to heating from disk is a weak 
function of a. Table [3] further indicates a ra- 
tio i?sub,disk/-Rsub,* 3 as a good estimate over 
a wide range of mass accretion rates and stellar 
mass. 

The value of a is related to the radial trans- 
port of matter in the disk, thus the mass accretion 
rate is a function of a. For Tabled the values of 
the sublimation radius is estimated using eq (|13p 
assuming the same mass accretion rate for differ- 
ent a values. These are representative values of 
sublimation radius for typical mass accretion rates 
found in massive star forming regions. Using these 



values we optimize the ratio of the self sublimation 
radius of the disk to that caused by heating from 
stellar luminosity is around three. With this as- 
sumption, we obtain a relation of a and M for a 
particular stellar mass. This is obtained by setting 
T — 1500 K in the equation and the radial 
distance as 3 x i?sub.*- 

This relation can be used to have constraints 
on the mass accretion rates for a particular stellar 
mass. The plot for mass accretion rates with a is 
shown in Fig. [2] for typical OB type stars. If one 
considers the value of a as fixed to 1, then from 
the Fig.[2]for the lOAf© star one gets the accretion 
rate of 4.2 x lO~^M0yr~^, which implies that for 
this order of mass accretion rate the dust in the 
disk will sublimate at a distance of 3 x i?sub.*- This 
clarifies that for the typical high accretion rates re- 
quired for formation of high mass stars, the heat- 
ing in the disk is very efficient to sublimate most 
of the dust grains in the midplane before the ra- 
diation could have any major effect on them. 

The implication of this result is profound, as 
it demonstrates that viscous heating of the disk 
is the dominant mechanism in the midplane for 
sublimation of the dust. The self-sublimation in 
the turbulent massive disk sublimates most of the 
dust grains well before the stellar radiation could 
affect them. Essentially, this implies that more 
matter (in form of gas) can reach closer to the 
central star. 

However, stellar radiation can have significant 
effects on the surface layer of the disk. The radial 
distance beyond which the stellar irradiation will 
dominate can be estimated by equating the flux 
from central star to the flux emitted from the disk. 
The flux from the disk is dependent on the effective 
surface temperature, which can be substantially 
lower than the midplane temperature depending 
on the optical depth, (see §4.3) 

4.3. Inner Gaseous disk structure 

In order to model the inner gaseous region of 
the disk, we consider the OPAL opacity tables ap- 
plica ble for higher temperatures ( Iglesias fc Roger j 
19961) . The various dynamical quantities in the 
disk obtained using the opacity from the table 
and their comparison with the various analytical 
models are shown in Fig. [3l The variation of the 
opacity, used for the present work, from the opac- 
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Table 3: Comparison of the sublimation radius due to heating in disk [equations & p^J and due t o 
heating from the star [equation (fTS]) ]. The temperature of the star and its radius are taken from Land (1992). 
(see §4.2 for details) 



Typical M [Mq yr"^] 


Spectral Type 


Stellar Mass (Mq) 




.[AU] 


-Rsub.disk 


[AU] 








Measure A 


Measure B 


a = 0.1 


a = 1 


10"^ 


B5 


5.9 


1.4 


1.0 


6.1 


3.5 


10-'' 


B2 
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ity tables with the midplane temperature in the 
disk is shown in Fig. [D 

The midplane temperature reaches very high 
values of the order of 10^ K at 2OR0 (See Fig. [3]). 
The radial profile for the temperature also shows a 
sudden break around 1 AU. This is related to the 
sudden rise in opacity around 3000 K as demon- 
strated in Fig. [4] Since the surface density and the 
mass density are inversely proportional to the tem- 
perature (see equations H]), ([3]), ([5]) & ([7])), their 
radial profiles show a sudden fall at that distance. 
Similar kind of sudden rise and fall can be seen 



also in the analytical results of iBell fc LinI (|1994| ) 
which are we show alongside in Fig [S] 

We also compare with analytica l radial profiles 
obtain ed using opacities as given bv lRuden fc Pollack 
(|l991l ). However these profiles just extent to the 
temperature range where the dust just sublimates 
and the opacity drops to a small value. These 
models fit very well in the outer disk region with 
the model described in the present work. 
The high temperature upto 10^ K would imply a 
high ionization fraction in the disk close to the 
massive star. Thus the ionized gas could well cou- 
ple to the large scale magnetic fields from ambient 
medium which is dragged by accretion and may 
give rise to collimated outflows. 

In Fig [21 the variation of the scale height is 
also shown with the radial distance. Similar to 
radial profiles of other dynamical quantities, we 
get jump also in the profile of the scale height. 
This jump is quite interesting when dealing with 
disks around massive stars. The sudden rise of 
disk height around 1 AU would help to shield outer 
disk regions from the radiation of central source. 
This affects the direction of radiation and leading 
to some sort of anisotropy in the radiation field. 



We also investigate the radial profiles of the dy- 
namical quantities for higher mass stars. These 
profiles are shown in Fig. [6l The midplane tem- 
perature profiles clearly show that as the central 
mass and accretion rate increase, the midplane 
temperature also increases and so the self sub- 
limation radius will move further out. For in- 
stance, with central mass as 23M0 and accretion 
rate of 4 x lO~'^M0yr-^, the dust sublimates at 
~ 30 AU whereas the dust sublimation radius is 
around 41 AU for a 37M0 star with disk accre- 
tion rate of 8.5 x lO~'^M0yr-^. These parameters 
are also chosen with the same argument as used 
for IOM0 star. Also, the temperature in the in- 
ner most region is much higher for massive young 
star with higher mass. The mass density and the 
surface density profiles also show the same trend 
of increment in the inner region with increase in 
mass of central object. 

On might question the thin disk approach to 
study the disk around massive young stars. How- 
ever, with this approach we find that the disk 
around high mass stars are very similar to that 
around cataclysmic variables. The spectral signa- 
tures of disk around these compact objects were 
very well explained by this thin disk model. The 
opacity values from the OPAL opacity table are 
used to study the bounda ry layer of the white 
dwarfs ( Collins et al.|[l998l ). They obtain a very 
high central temperature near to the white dwarf 
but on applying the optical depth consistently 
from the table, they get the same order of magni- 
tude of the effective temperature 



gy4 
"37 



as it is on the surface of a typical white dwarf. 
In our case, the midplane temperature profile 
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shows that the region very near 2OR0) is very 
hot gas of 10^ K, but to compare with the stellar 
surface effective temperature (Toft = 22000 K for 
M = IOM0) one has to take into account the finite 
optical depth and then the surface temperature of 
the disk, plotted in the same figure matches the 
surface temperature of the star at the very inner- 
most regions. Thus the effective temperature of 
the disk can be considerably lower than the mid- 
plane temperature. The flux from the disk and the 
star are equated in order to compute the radial dis- 
tance where the stellar irradiation can affect the 
disk surface, 



47rr^ 



off 



(16) 



where L^, is the luminosity of the star and cr is 
the Stefan-Boltzmann constant . We find for a 
IOM0, the stellar irradiation will have a effect on 
the surface layers around ^ 30 AU for a = 1. If 
the a value is lowered, the stellar irradiation can 
have affect on the surface at even smaller radial 
distance of around 10 AU. However this effect will 
not alter the dust sublimation radius in the disk. 

4.4. Growth of Massive Stellar Embryo 

One of the major problems in the formation of 
massive stars is the large UV radiation from the 
protostars that exert pressure on the matter in- 
hibiting the infall on them. This is because as the 
central mass increases also the central luminosity 
increases, thus exerting a large radiation pressure 
which is usually thought to halt the matter falling 
on the star. In order to get away with th e large 
radiant flux and to form more massive stars, lYorkd 
lists some of the favorable conditions such 
as i) Reduction of Kes ii) Reduction of effective 
Luminosity and Hi) Increasing the gravitational 
acceleration. 

Reduction of k^s and effective Luminosity 

In earlier approaches considering the conditions 
for the formation of massive stars, the matter was 
thought to be accumulated by spheric al accretion 
(|Kahnlll974l : [Wolfire fc Cassinellilll987l ). Using de- 
tailed dust models and complex grain size distri- 
bution, these models were able to put a limit on 
the maximum stellar mass of 20Mo. In general, 
the dominance of the radiation pressure over grav- 
ity of the spherical infall of matter puts a strong 



constraint on the value of the opacity. In this case, 
the necessary condition to form a massive star is 



(17) 



fe.g. lYorkell2004l : IZinnecker fc Yorkell2007n . which 
implies an upper limit to the effective opacity. 



Koff < ISOcm^g ^ 



M 

lOMr. 



L 



lOOOLp 



(18) 



In the disk scenario, there would be additional 
contribution to the forces mentioned in equation 
(fT7| . such as disk gas and ram pressure. 

In the present work, we demonstrate that most 
of dust, due to the self sublimation of accretion 
disk, is destroyed already at distances larger that 
~ 10 AU from the central source which helps to 
reduce the Keff- It is evident that the maximum 
value of the opacity even in the innermost region 
of the disk is less than 5 cm^g^^ (see Fig. |4]), 
which is much smaller than the upper limit given 
by equation (|18|) required for the formation of a 
typical massive star. This suggests that the radia- 
tion pressure is no longer inhibiting the accretion 
process onto to the central star and thus allowing 
more massive stars to form. This effect is aid- 
ing the matter to flow closer to the central object 
but does not ensure that matter will be accreted 
on the central star. Thus reduction in opacity is 
an essential requirement but it does not guarantee 
growth of mass as the isotropic radiation from the 
star may still stop the dust by exerting the pres- 
sure which is dependent on the luminosity of the 
central object. 

A sketch of our model calculations assuming a 
IOM0 central star is shown in Fig. [71 The scal- 
ing of the figure is done using a fiducial jet radius. 
The main components in the figure are the dusty 
disk which is feeded with matter by the core, the 
inner gaseous disk and the large scale bipolar out- 
fiows. The region near to the dust sublimation 
radius shows the presence of dust and gas and 
large variation of opacity is seen in this region. 
Very close to the source there is large fiux of UV 
and visible radiation which is blocked from the 
midplane of outer dusty disk by the high optical 
depths (r ^ 1) of the inner gaseous species. The 
figure also shows the possibility of radiation flux 
escaping more in the direction perpendicular to 
the plane. This "flashlight effect" was introduced 
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bv lYorke fc Bodenheimed (tl999h in order to have 
anisotropy in the radiation and allow matter to 
accrete. 

4-4-^- Overcoming Radiation Pressure 

The reduction of effective opacity and luminos- 
ity ensures the matter to come closer to the cen- 
tral star. In the region very near to the star mat- 
ter is subjected to different pressure sources, all of 
them play an important role in understanding the 
dynamics of disk accretion under the influence of 
stellar radiation pressure. We therefore compare 
the contribution of different pressure sources, con- 
sidering typically a IOM0 star. 

The radiation pressure from the star is a func- 
tion of the radial distance from the surface of the 
star, 



Ancr"^ 



= 5.9 X lO^ergcm ^ ' 



cff 



V 22000K 



57.1 



The ram pressure from the disk counteracts the 
radiation force from star and will allow the mat- 
ter to accrete. We estimate the ram pressure 
-Pram = pv^ using our consistent dust and gas 
opacity model (see §4.1 and §4.3), 



P,^^{r) = 2.08 X lO^ergcm" 



M 



P 



2.32 X 10-9gcm- 



M 



T 



lO-^Moyr-i / VlOM©/ V5.03x lO^K 



lOR* 



Figure [8] shows the radial profiles of both pres- 
sure sources. The figure clearly depicts that the 
radiation pressure from the star is much below the 
ram pressure in the disk which will aid the matter 
to overcome radiation pressure and accrete on the 
central star. In the very inner region of the disk, 
the disk radiation pressure becomes comparable to 
the gas pressure in the disk which may result in 
thermal instability. (Figure [5]) 

5. Stability of Disks 

The disks around massive stars can be unsta- 
ble due to axisymmetric gravitational instability. 



They can also be unstable due to fragmentation 
due to rapid cooling of the disk as compared to the 
dynamical time sca le. ( Gammid 2001 : Rice et al 



2003t iRafikovllioOSI ). We will also see that these 



disks can be also unstable thermally in the very 
inner region, close to the massive star. 

5.1. Gravitational Instability 

The criterion for the gravitational instability of 
a disk is given by the Toomre parameter 



Q 



27rGE 



with as the Kcplcrian angul ar velocity and E as 
the surface density of the disk ("Toomre" 1964) . For 
Q > 1 the disk is stable (no fragmentation), while 
Q <1 leads to instability in the disk. Physically, 
the Toomre parameter can be considered as the ra- 
tio of centrifugal force along the radial direction to 
the gravitational force acting in the direction per- 
pendicular to the radial motion. Thus, if a local 
accumulation of mass is moving in a certain orbit, 
then this would lead to slowing down of orbital 
motion and also more gravitational force acting 
downward which implies that the value Q will de- 
crease and eventually when Q < \ the downward 
gravity force wins and leads to instability gener- 
ating overdense regions. 



It has be en known from simulations (jKrumholz et al 

\2mi. '2009') and observatio ns in case of G 192.16- 
-1 3.82 (Shepherd et al.l l200ll . considering large ob- 
servational errors) that disks in massive stars are 
unstable as they may have Q value low and some- 
times < 1 for some radial extent. 
(20) rpj_^g Toomre parameter decreases with the ra- 
dial distance, and at some radial distance and 
the disk becomes gravitational instable. Figure |4] 
shows the radial behavior of the Toomre param- 
eter in case of typical disk parameters (M = 
4.2 X 10"'' Mq yr-i , stellar mass M = IOMq 
and a ^ 1). It is clear that such a disk is unstable 
as g < 1 after - 100 AU. 



tThe Toomre criteria can be written as a product of two 
ratios - -^tfeermai E^„t <; ^ jf onlv one of the ratios is 
less than unity, that will not guarantee that fragmentation 
will occur, (see §5.2) 
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5.2. Fragmentation of Disks 

Disks which are stable for axisymmetric grav- 
itational instability may not be necessarily sta- 
ble for fragmentation which leads to formation of 
bound objects. Gravitational instability sets an 
upper limit of the sound speed where as the analy- 
sis for fragmentation and cooling time set s a lower 
limit on the speed of sound in the disk! Rafikovl 
20051) . 

Nu merical simulations ( Gammiell2001 [ Rice et al 



20031 ) and analytical solutions (|Rafikovll2005f ) sug- 
gest that for a disk to avoid fragmentation, the 
cooling time scale should be larger than the dy- 
namical timescale {tdyn ~ fi^^). The thresh- 
old relation of the cooling time can be obtained, 
^tcooi < Ci where the factor C « 3, as obtained 
from simulations usin g constant tcooi (|Gammie 
20011: iRice et al.|[2003h . For the present purpose, 
we apply cooling time as given bv lRafikovl feoOSl 



^cool 



^cl f{r) 



1- 



1 2crT4 



(21) 



where /(r) — t + - describes the efficiency of 
cooling and depends on the optical depth r. In 
the above equation ([2T|) . 7 is the adiabatic index 
and f actor 7 — 1 considered to be of the order of 



unitv (|Rafikovll2005D . 

The variation of the cooling time, as given by 
eq. ([2T|l . with radial distance is shown in Fig. |4] 
for different values of a. These curves are made 
applying the Bell & Lin opacity power laws. The 
solid line in the figure for cooling time determines 
the threshold for fragmentation to set in. The re- 
lation suggests that a value a ^ 1 would lead to a 
disk which may fragment completely. However, 
the lower values of a ~ 0.1,0.01 would enable 
a stable disk against fragmentation. The value 
a ~ 0.1 is consistent with that obtai ned by ob- 



serva tional modelling of ionized disks (jKing et al 
20071 ). 



5.3. Thermal Instability 

Thin accretion disks around black holes which 
are supported by radiation pressure (region A, 
See §2) were found to be thermally instable 
(IShakura fc SunvaevI [l976l) . The physical reason 
for such kind of instability is inefficient cooling in 
the disk as compared to the viscous heating. This 
leads to overheating which causes expansion which 



in turn overheats the disk eventually leading to a 
thermal runaway. There were ideas of application 
of slim disk model to this region, which is a stable 
branch in the S-shaped M — E curve, because such 
a model prevents the radiation to escape from the 
disk and allowing the flux to be advected alon g 
with the flow of matter |Abramowicz et allllQSSl ). 

In the modelling of the disk around massive 
star, we see that the radiation pressure becomes 
comparable to the gas pressure in the very inner- 
most region of the disk < AR^, for a lOM© star. 
Here the assumption that Pgas > ^'rad no longer 
holds, and we stop our iterations at this radii. 
However, the application of concept of slim disk 
in protostellar disks for efficient cooling is not re- 
ally probable. This is because of the fact that 
advected flux should get rid in some manner, it is 
viable in disks around black holes as the flux can 
be advected into the black holes, however in case 
of disks around stars this flux will heat up the star , 
which may be not physical. iTurner et al.l (|2007h 
applied the concept of photon bubble instabilities 
to the young massive stellar environments, which 
could be an efficient way of cooling the innermost 
region of massive star forming region as well. 

5.4. How big is al 

The value of a parameterizes the viscosity in 
the disk. In the present case, we see that higher 
value of a ~ 1 ensures that the disk within ~ 
lOOAU is stable to axisymmetric gravitational in- 
stability, though for this high value of a , the disk 
completely would fragment. 

For a ~ 0.1 the disk would become stable to 
fragmentation. In this case, the part of the disk 
stable to axisymmetric gravitational instability re- 
duces to < 80 - 90AU 

For a — 0.01, the disk would again be stable 
to fragmentation, but now most part of the disk 
is subjected gravitational instability. The effect of 
lowering the value a on the Toomre parameter can 
be seen in Fig. |4] 

The value of a is usually chosen as a param- 
eter in disk modelling for different purposes. In 
case of low mass stars, the range of a values 
accepted is around 0.001-0.01, which is believed 
to be produced by mag neto-rotational instability 
(|Balbus fc Hawlevlll998l ). In case of massive stars, 
the value of a can be higher by one order of mag- 
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nitude and the physics involved to generate this 
high a value could be different from that seen in 
low stellar mass case, (for e.g. Gravitational in- 
stability) 

A high a would imply that matter flows ra- 
dially with high speed near to local sound speed 
and there could be a possibility of accretion shocks 
produced to make the radial speed of matter sub- 
sonic. These accretion shocks have been proposed 
for the close circumstellar environments of high 
mass p rotostar such as CRL 2136 w ith the H2O 
maser (jMenten fc van der Takll2004 ). 

The other way of treating viscosity in the disk 
would be by assuming a radial profile, a(r), such 
that its value may be high (~ 1) in the outer re- 
gion and lower(~ 0.1) in the inner region of the 
disk. Such a radial profile would then lead to frag- 
mentation resulting in formation of bound objects 
influencing the structure in the outer part of the 
disk around most massive star formed such that 
each bound object may form itself a low m ass or 
even high mass star ( Krumholz et al. 2009f ). How- 
ever, the inner ionized disk can be still stable and 
can launch large scale bipolar outflows. 

In Fig. [51 variation of the dynamical quanti- 
ties obtained from the present model for a lOM© 
star with two different values of a is shown. Apart 
from changing the a value we also change the value 
of mass accretion rate according to the figure [2] 
There is no considerable change in midplane tem- 
perature profile, which is due to the fact that by 
changing two parameters also bring change in the 
opacity and the cumulative effect cancels the vari- 
ation in temperature [§. However, this does not 
happen for mass density as the dependence on the 
above variations is different. Since the height de- 
pends only on the temperature, the height profile 
also does not show any variation. 

6. Implications for Observations 

Early stages of massive star formation are dif- 
ficult to observe, as the whole process is enclosed 
in an envelope of dust. Also these high mass star 
forming regions are at a distance of few kilo par- 
sees which makes the observation more difficult 
due to limited resolution of the present day tele- 



scopes. 

One of the essential findings of the present work 
is the high temperature of the order of ~ 10^ K in 
the inner disk. At such a high temperature, the 
opacity is mostly dominated by electron scattering 
and there is a pool of highly energetic electrons 
present close to the massive star. These electrons 
might be responsible for soft X-ray emissions from 
massive young stars. 

There have been observations of X-rays in some 
of th e massive young so urces in the W3 com- 



plex dHofner et all [iooi l and also in GGD 27 



(jPravdo et al.ll2009f ). In the case of GGD 27, ob- 
served hard X-rays emissions have energies of the 
order of 2-10 KeV. The temperatures obtained 
from present work indicate emissions of soft X- 
rays, however the soft X-rays emissions have a 
large extinction in massive star forming regions. 
There might be possibility that these less energetic 
electrons are accelerated by magnetic fields giv- 
ing rise to hard X-ray emission and also the hard 
X-ray emission may come from shocks produced 
when the accreted matter falls on the surface of 



star (|Montmerldl2007n . 

Amidst all the difficulties, observers were suc- 
cessful to resolve a fattened structure having 
length scales of ~ 1000 AU. The best candidate of 
showing presence of a rotating Keplerian disk is 
IRAS 21026 -I- 4104 (|Cesaroni et al.ll2005h which 
is predicted to be early B star (Mass = 7-12 M©). 
They estimate the radial profile of mean kinetic 
temperature for the Keplerian disk and mass den- 
sity profile using different lines. They obtain the 
best fit for temperature profile as T{R) oc R^"^ 
where q is 0.57-0.75 and that of mass density 
p{R) = i for HCO+ (1-0) line. The val- 
ues of the powers obtained are very close to the 
best fit o f these profiles from the present model 
(Fig. H]). ISchrever et all (|2006h also were able to 
resolve a rotating structure in an embedded 8-10 
Mq massive star AFGL 490. They also obtain a 
Keplerian profile for the rotational velocity and 
the best fit radial power law for the surface den- 
sity was E cx R~^'^, which is close to the value 
obtained from the present model (E oc R~^'^). 
In case of low mass stars, usually the density in 
the disk is fitted using the gaussian as given by 



^Note that S is inversely proportional to a and decreasing 
the a will increase density and thus opacity 
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()Whitnev et al.ll2003l : IWolf et al.ll2008h 

/i?.V J-l ' - 

Pdisk = I ~ ) ^^P]^ 



(22) 



where the scale height has the radial dependence 
of the form 



(23) 



The best fit obtained for the cas e of disks in "But - 
terfiy" star was a=2.35, b=1.28 (|Wolf et al.ll2008f) . 
These results are interesting as similar radial fits 
to the elongated continuum are obtained even for 
the recent observational results of young high mass 
source IRAS 18151-1208 (Fallscheer private com- 
munication). Also similar fits are obtained from 
the present model. 

7. Conclusion 



In this paper we have presented - for the first 
time - global models for accretion disks around 
massive young stars. 

We have solved the thin disk equations tak- 
ing into account the proper opacities for dust and 
gas. In particular we consider Rosseland mean 
dust opacities given bv iRuden fc Pollack ( 1991 ) 
and the gas opacities from t he OPAL opacity ta- 
bles ( Iglesias fc Rogerslll996 ). 

This enables us to provide the dynamical quan- 
tities of disk accretion from the very inner part 
at radii of 0.1 AU to the outer region of the disk 
at about 100 AU. At the same time this provides 
a theoretical link between outer disk which is ac- 
cessible in principle by observations and the inner 
disk which is not yet possible to resolve observa- 
tionally. 

Our main results can be summarized as follows. 



3. We estimate the stability of these disk by the 
Toomre criterion and find that our thin disks 
around e.g. a 10 Mq star becomes gravita- 
tionally unstable beyond 100 AU. We also 
see for a values close to 1, the disk would 
fragment completely whereas for a ~ 0.1 the 
disk could remain stable to fragmentation. 
We also discuss the effect on the dynamics 
and stability of the disk model with the vari- 
ation of a. 

4. For the given disk and stellar parameters and 
disk opacities we find that the stellar radi- 
ation pressure is negligible against the disk 
ram pressure and gas pressure and therefore 
cannot hinder accretion towards a massive 
young star. 

Considering the high disk temperatures and the 
rather large disk viscosity parameter, disk around 
massive young stars seem to be intrinsically differ- 
ent from the low mass equivalents in particular to 
the different form of disk opacity 

We also find that for higher mass stars with 
high accretion rates the dust sublimates further 
away and obtain much higher temperature for the 
same. The presence of the optically thick gaseous 
component and anisotropy of radiation prevents 
the matter to be unaffected by the radiation from 
the star. 



1. For typical stellar masses and accretion rates 
we find very high midplane temperatures of 
the order of 10^ K for radii less than ^ O.lAU 

2. Due to the high disk temperatures the dust 
sublimates already at distances which are 
about a factor three larger than caused 
by the stellar irradiation. This disk self- 
sublimation lowers the disk opacities consid- 
erably and allows for disk accretion in the 
stellar radiation field. 
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Fig. 1. — Radial profiles of midplane temperature (T), mass density (p), surface density (S) and the scale 
height ratio {h/r) for the outer dust dominated disk. The dashed l ine re presents the tabulated values (see 
Table [2|), the solid line is obtained by using the iRuden k, Pollackl ( 199l[ ) opacity power laws. These plots 
are for typical M = 4.2 x 10"'' yr~^ , stellar mass M = IOM0 and a ^ 1. The radial index of the best 
fitted linear plot are mentioned in each subplot. 
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Fig. 2. — This plot shows the variation of mass accretion rate with the viscosity parameter so that the 
dust subhmation temperature due to viscous heating is around three times that of heating from star. The 
various hues are for different spectral type of stars (luminosity values obtained from iLane..l992) . These plots 
indicate that mass accretion rate is a weak function of a. 
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scale h eight of the disk exte nded to the 
I2OO5I) and O PAL 



Fig. 3. — Radial profiles of temperature (T), mass density (p), 

i nner region. The solid line represents the values from opacity table (iFerguson et al. 
(jlglesias fc Rogers , the dashed line is using the opacity power laws bv' Ruden fcPoUackl (|l99lh and 

the dot-dashed line are the values for Bell and Lin opacity power laws. The dotted line shown in temperature 
profile represents the T^e profile with radius whereas the other lines are for the midplane temperature. These 
plots are for typical M = 4.2 x lO""* M© yr'^ stellar mass M — lOM© and a ^ 1 
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Fig. 4. — The opacity variation with temperature [top left panel), and the variation of Toomrc parame- 
ter Q with the radia l distance (t op right panel). The so lid line represents the values from opacity table 
((Ferguson et al."2005') and OPAL (jlglesias fc Roger the dashed line is using the opacity power laws 

bv .Ruden fc P ollack (1991) and the dot-dashed line are the values for Bell and Lin opacity power laws. The 
bottom panels shows how the variation of a values affect the cooling time and Toomre parameter, a — 0.01 
is denoted by dotted line, where as a = 0.1 is dashed line and a = 1.0 is by the dot-dashed line. Solid line 
in the curve represents the threshold condition for the fragmentation to set in. 
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Fig. 5. — Dynamical quantities for two different a values. The solid line represents plot with a = 0.1 and 
M = lO"** where as dashed line is for a = 1 and M = 4.2 x lO^** 



18 



10« 



1^ 



T — I I I I Mil 1 — I I I Mill 1 — I I I I INI 1 — TT 




2 Q3 I I I I I I I I L_L 

0.1 1 10 100 



lO-* r 



n 

e 

o 
\ 

60 




10- '3 



'III iii'l ' I l_L 



0.1 1 10 100 

R (AU) 



-| — I I I I INI 1 — I I I Mill 1 — I I I I INI 1 — TT 




J Ql I I I I I I I I L_L 

0.1 1 10 100 

R (AU) 



M = 10 M„ 



M = 23 M„ 



- M = 37 M« 



Fig. 6. — Radial profiles of temperature (T), mass density (p), surface density (S) of the disk extended to 
the inner region. The solid line corresponds to M = lOM©, the dotted line for M = 23M0 and the dashed 
line is for central stellar mass M = 37Mq. 
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RADIATION CARRIED AWAY VIA CAVITY 

* 
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Fig. 7. — Pictorial representation of the inner region of the massive star. 
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Fig. 8. — The above figure compares the pressure due to different sources in the circumstellar environment 
for a 10 Mq star 
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